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We present a precise measurement of the ratio Rk — T(K — > ev{y))/T{K — > juv(y)) performed 
with the KLOE detector. The results are based on data collected at the Frascati e + e~ collider 
DA"NE for an integrated luminosity of 2.2 ftr 1 . We find R K = (2.493 ±0.025 sta t±0.019 syst ) x 
10~ 5 , in agreement with the Standard Model expectation. This result is used to improve con- 
straints on parameters of the Minimal Supersymmetric Standard Model with lepton flavor viola- 
tion. 
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The decay K —>■ e ± v is strongly suppressed, ~fewx 10~ 5 , because of conservation of angular 
momentum and the vector structure of the charged weak current. It therefore offers the possibility 
of detecting minute contributions from physics beyond the Standard Model (SM). This is partic- 
ularly true of the ratio R K = T(K — ► ev)/T{K — ► }iv) which, in the SM, is calculable without 
hadronic uncertainties [jl|, g]. Physics beyond the SM, for example multi-Higgs effects inducing 
an effective pseudo-scalar interaction, can change the value of Rk- It has been shown in Ref. [||] 
that deviations of Rk of up to a few percent are possible in minimal supersymmetric extensions of 
the SM (MSSM) with non vanishing e-x scalar lepton mixing. To obtain accurate predictions, the 
radiative process K — > evy (K e 2y) must be included. In K e %y, photons can be produced via internal- 
bremsstrahlung (IB) or direct-emission (DE), the latter being dependent on the hadronic structure. 
Interference among the two processes is negligible Pp. The DE contribution to the total width is 
approximately equal to that of IB [Q] but is presently known with a 15% fractional accuracy [§]. 

Rk is defined to be inclusive of IB, ignoring however DE contributions. A recent calculation 
[g|, which includes order e 2 p 4 corrections in chiral perturbation theory C^PT), gives: 

R K = (2.477 ±0.001) x 1(T 5 . (1) 

Rk is not directly measurable, since IB cannot be distinguished from DE on an event-by-event 
basis. Therefore, in order to compare data with the SM prediction at the percent level or better, one 
has to be careful with the DE part. We define the rate Rio as: 

R w = r(K^ev(y), E 7 < 10 MeV) /V (K ^ nv). (2) 

Evaluating the IB spectrum to &(a em ) with resummation of leading logarithms, Rio includes 
93.57 ±0.07% of the IB, 

R w = R K x (0.9357 ± 0.0007) . (3) 

The DE contribution in this range is expected to be negligible. Rio is measured without photon 
detection. Some small contribution of DE is present in the selected sample. In particular, DE 
decays have some overlap with the IB emission at high p e . We have also measured [0] the dif- 
ferential width dR r /dEy for E r > 10 MeV and p e > 200 MeV requiring photon detection, both to 
test #PT predictions for the DE terms and to reduce possible systematic uncertainties on the Rio 
measurement. 

1. Selection of lep tonic kaon decays 

^ decays are signaled by the observation of two tracks with the following conditions. One 
track must originate at the interaction point (IP) and have momentum in the interval {70, 130} 
MeV, consistent with being a kaon from <p -decay. The second track must originate at the end 
of the previous track and have momentum larger than that of the kaon, with the same charge. 
The second track is taken as a decay product of the kaon. The point of closest approach of the 
two tracks is taken as the kaon decay point D and must satisfy 40< ro <150 cm, |zd| <80 cm. 
The geometrical acceptance with these conditions is ~56%, while the decay point reconstruction 
efficiency is ~51%. From the measured kaon and decay particle momenta, and pa, we compute 
the squared mass mj of the lepton for the decay K — > iv assuming zero missing mass: 

mf=(£*-|P*-Pd|) 2 -Pd- (1.1) 
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The distribution of m 2 is shown in Fig. [I] left panel, upper curve, from MC simulation. The muon 
peak is quite evident, higher masses corresponding to non leptonic and semileptonic decays. No 
signal of the K — > ev decay is visible. The very large background around zero mass is the 
tail of the K — > /xv (^2) peak, due to poor measurements of pk, Pa or the decay angle, and- The 
expected signal from K e %y is also shown in Fig. [j] left, lower curves, separately for E y >10 and 
<10 MeV. The expected number of K e 2 decays in the sample is ~30,000. A background rejection 
of at least 1000 is necessary, to obtain a 1% precision measurement of Y{K e 2), with an efficiency 
of ~30%. 




m\ (MeV 2 ) m 2 , (MeV 2 ) 

Figure 1: Left: MC distribution of mj, solid line. The contribution of K&y with E r < 10 MeV (> 10 MeV) 
is shown by the dashed (dotted) lines. Right: mj spectrum before (dashed) and after (solid) quality cuts for 
MC K^2 (upper plots)and K e 2 with Ey < 10 MeV (lower plots). Black dots represent data after quality cuts. 

The kinematics of the two-body decay <p — * K + K provides an additional measurement of pk 
with ~ 1 MeV resolution, comparable with that from track reconstruction. We require the two px 
determinations to agree within 5 MeV. Further cuts are applied to the daughter track. Resolution of 
track parameters is improved by rejecting badly reconstructed tracks, i.e., with % 2 (track fit)/ndf > 
7.5. Finally, using the expected errors on pk and p& from tracking, we compute event by event 
the error on rrq, Srrq. The distribution of 8m q depends slightly on the opening angle Oga, which 
in turn has different distribution for K e % and K^%. Events with large value of 8mj are rejected: 
Srrq < 5 max , with cV ax defined as a function of OCxd, to equalize the losses due to this cut for K e i 
and Ku2. The effect of quality cuts on m 2 resolution is shown in Fig. Ill right. The background in 
the K e 2 signal region is effectively reduced by more than one order of magnitude with an efficiency 
of ~70% for both K e2 and K^. 

Information from the EMC is used to improve background rejection. To this purpose, we ex- 
trapolate the secondary track to the EMC surface and associate it to a nearby EMC cluster. This 
requirement produces a signal loss of about 8%. Energy distribution and position along the shower 
axis of all cells associated to the cluster allow for ej\l particle identification. For electrons, the 
cluster energy E c \ is a measurement of the particle momentum p<], so that E c i/pd peaks around 1, 
while for muons E c \/ pj is on average smaller than 1. Moreover, electron clusters can also be dis- 
tinguished from jit (or %) clusters by exploiting the granularity of the EMC. All useful information 
about shower profile and total energy deposition are combined with a 12-25-20-1 structure neural 
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network trained on Kl — ► 7l£v and K^i data, taking into account variations of the EMC response 
with momentum and impact angle on the calorimeter. The distribution of the neural network output, 
AW, for a sample of Kl — > nev events is shown in Fig. || left, for data and MC. Additional separa- 
tion has been obtained using time of flight information. The data distribution of AW as function of 




Figure 2: Left: Neural-network output, AW, for electrons of a Kl — > Kev sample from data (black) and MC 
(red). Right: Data density in the AW, w» plane. 

m\ is shown in Fig. fright. A clear K — » ev signal can be seen at ra 2 ~ and WW ~ 1. 

Some 32% of the events with a K decay in the fiducial volume, have a reconstructed kink 
matching the required quality criteria and an EMC cluster associated to the lepton track; this holds 
for both K e 2 and Ku2- In the selected sample, the contamination from K decays other than K(2 is 
negligible, as evaluated from MC. R\q, Eq. ^[ is obtained without requiring the presence of the 
radiated photon. The number of K — > ev(y), is determined with a binned likelihood fit to the two- 
dimensional WW vs ra 2 distribution. Distribution shapes for signal and Ku% background are taken 
from MC; the normalization factors for the two components are the only fit parameters. The fit has 
been performed in the region -3700 < mj < 6100 MeV 2 and AW > 0.86. The fit region accepts 
~ 90% of K — > <?v(y) events with E y < 10 MeV, as evaluated from MC. A small fraction of fitted 
K — > ev(y) events have Ey > 10 MeV: the value of this "contamination", /de> is fixed in the fit 
to the expectation from simulation, /de = 10.2%. A systematic error related to this assumption is 
discussed in Sect. ^. 

We count 7064±102 K+ -> e+v(y) events and 6750±101 K -» e~v(y), 89.8% of which 
have Ey < 10 MeV. The signal-to-background correlation is ~ 20% and the # 2 /ndf is 113/112 for 
K+ and 140/112 for K . Fig. | shows the sum of fit results for K + and K projected onto the m 2 
axis in a signal (WW > 0.98) and a background (WW < 0.98) region. 

The number of K^ 2 events is obtained from a fit to the mj distribution. The fraction of back- 
ground events under the muon peak is estimated from MC to be less than one per mil. We count 
2.878 x 10 s (2.742 x 10 8 ) K + -> Ai + v(y) (K~ -> ^v(y)) events. The difference between K + and 
K counts is due to K nuclear interactions in the material traversed. 

2. Efficiency and systematic errors evaluation 

The ratios of K e 2 to efficiencies are evaluated with MC and corrected for possible differ- 
ences between data and MC, using control samples. We evaluate data-MC corrections separately 
for each of the following analysis steps: decay point reconstruction (kink), quality cuts, cluster- 
charged particle association. For each step, the correction is defined as the ratio of data and MC 
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Figure 3: Fit projections onto the mj axis for AW > 0.98 (left) and AW < 0.98 (right), for data (black dots), 
MC fit (solid line), and background (dotted line). The contribution from K e 2 events with Ey > 10 MeV 
is visible in the left panel (dashed line). 

efficiencies measured on the control sample, each folded with the proper kinematic spectrum of 
K e 2 (or K^z) events. 

Decay point reconstruction efficiencies are evaluated using pure samples of an(1 these 
are tagged by the identification of the two-body decay, or K — ► 7171° (K^), of the other kaon 
and selected with tagging and EMC information only, without using tracking. The corrections to 
MC efficiencies range between 0.90 and 0.99 depending on the decay point position and on the 
decay angle. The simulation is less accurate in case of overlap between lepton and kaon tracks, and 
with decays occurring close to the inner border of the fiducial volume. 

Samples of K L (e3), K L (n3), and K^ 2 decays with a purity of 99.5%, 95.4%, and 100.0% 
respectively, are used to evaluate lepton cluster efficiencies. These samples are selected using 
tagging and DC information only, without using calorimeter, see Refs. [§J. The efficiency has been 
evaluated as a function of the particle momentum separately for barrel and endcap. The correction 
to MC efficiencies ranges between 0.98 and 1.01 depending on the momentum and on the point of 
impact on the calorimeter. The trigger efficiency has been evaluated solely from data. 

The absolute values of all of the systematic uncertainties on Rio are listed in Table [I]. All of 
the sources of systematic error are discussed below. 



5(/?io)xl0 5 


Statistical error 


0.024 


Systematic error 




Counting: fit 


0.007 


DE 


0.005 


Efficiency: kink 


0.014 


trigger 


0.009 


e,/J. cluster 


0.005 


Total systematic error 


0.019 



Table 1: Summary of statistical and systematic uncertainties on the measurements of Riq. 

To minimize possible biases on K e % event counting due to the limited knowledge of the mo- 
mentum resolution, we used data to carefully tune the MC response on the tails of the rrn 
distribution. This has been performed in sidebands of the AW variable, to avoid bias due to the 
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presence of K e 2 signal. Similarly, for the NN distribution the EMC response in the MC has been 
tuned at the level of single cell, using data control samples. Residual differences between data 
and MC K e % and NN shapes have been corrected by using the same control samples. Finally, 
to evaluate the systematic error associated with these procedures, we studied the variation of the 
results with different choices of fit range, corresponding to a change of overall purity from ~ 75% 
to ~ 10%, for K — > ev{y) with Ey < 10 MeV. A systematic uncertainty of ~ 0.3% is derived by 
scaling the uncorrected errors so that the reduced % 2 value equals unity (see also Table [[]). 

K e 2 event counting is also affected by the uncertainty on /de> the fraction of K e 2 events in the fit 
region which are due to DE process. This error has been evaluated by repeating the measurement of 
R\q with values of /de varied within its uncertainty, which is ~ 4% according to our measurement 
of the K e 2y differential spectrum Since the distributions for K&y with Ey < 10 MeV and 
with Ey > 10 MeV overlap only partially, the associated fractional variation on /?io is reduced: the 
final error due to DE uncertainty is 0.2% (Table |]). 

Different contributions to the systematic uncertainty on £ e 2/ e n2 we listed in Table [[} These 
errors are dominated by the statistics of the control samples used to correct the MC evaluations. 
In addition, we studied the variation of each correction with modified control-sample selection 
criteria. We found neglible contributions in all cases but for the kink and quality cuts corrections, 
for which the bias due to the control-sample selection and the statistics contribute at the same level. 

The total systematic error is ~ 0.8%, to be compared with statistical accuracy at the level of 
~ 1%. 



3. Rk and lepton-flavor violation 

The number of K — > ev(y) events with Ey < 10 MeV, the number of K — > juv(y) events, the 
ratio of K e 2 to efficiencies and the measurement of R w are given in Table ^ for K + , K and both 
charges combined. K + and K~ results are consistent within the statistical error. The systematic 
uncertainty is common to both charges. 





N(K e2 ) 


N(^ 2 ) 




Rio x 10 5 


K+ 

K 

K± 


6348 ±92 ±23 
6064 ±91 ±22 
12412 ±129 ±45 


2.878 x 10 8 
2.742 x 10 8 
5.620 x 10 8 


0.944 ±0.003 ±0.007 
0.949 ± 0.002 ±0.007 
0.947 ± 0.002 ±0.007 


(2.336 ±0.033 ±0.019) 
(2.330 ±0.035 ±0.019) 
(2.333 ±0.024 ±0.019) 



Table 2: Number of K e 2 and events, efficiency ratios and results for R\q for K + , K 
combined; first error is statistical, second one is systematic. 



and both charges 



To compare the R\o measurement with the inclusive Rk prediction from SM, we take into 
account the acceptance of the 10 MeV cut for IB, Eq. |. We obtain: R K = (2.493 ± 0.025 stat ± 
0.019 syst ) x 10~ 5 , in agreement with SM prediction of Eq. []J In the framework of MSSM with 
lepton-flavor violating (LFV) couplings, Rk can be used to set constraints in the space of relevant 
parameters, using the following expression [||] : 



Rk 



■ R K M x 



i+l -r 

m H 



^c) U31I 



tan 6 ^ 



(3.1) 
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where M# is the charged-Higgs mass, A^ is the effective e-x coupling constant depending on 
MSSM parameters, and tan/3 is the ratio of the two Higgs superfields vacuum expectation values. 
The regions excluded at 95% C.L. in the plane M//-tanj3 are shown in Fig. [| for different values 
of the effective LFV coupling A^ 1 . 




200 400 600 800 
charged Higgs mass (GeV ) 



Figure 4: Excluded regions at 95% C.L. in the plane M//-tan/3 for A 3 R l = 1CT 4 ,5 x 10~ 3 , 1CT 3 . 



4. Conclusions 

We have performed a comprehensive study of the process K e 2y. We have measured the ratio of 
K e 2y and Kui widths for photon energies smaller than 10 MeV, without photon detection require- 
ment. We find: R w = (2.333±0.024 st at±0.019 st at) X 1(T 5 . From this result we derive the inclusive 
ratio R K to be compared with the SM prediction: R K = (2.493 ±0.025 stat ±0.019 syst ) x 10~ 5 , in 
excellent agreement with the SM prediction Rk = (2.477 ±0.001) x 10~ 5 . Our result improves the 
accuracy with which Rk is known by a factor of 5 with respect to the present world average and 
allows severe constraints to be set on new physics contributions in the MSSM with lepton flavor 
violating couplings. 
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